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Abstract 

We present a geometric description of lepton flavor mixing and CP violation 
in matter by using the language of leptonic unitarity triangles. The exact 
analytical relations for both sides and inner angles are established between 
every unitarity triangle in vacuum and its effective counterpart in matter. The 
typical shape evolution of six triangles with the terrestrial matter density is 
illustrated for a realistic long-baseline neutrino oscillation experiment. 
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I. INTRODUCTION 



Recent solar [1], atmospheric [2], reactor (KamLAND [3] and CHOOZ [4]) and accelerator 
(K2K [5]) neutrino oscillation experiments have provided us with very convincing evidence 
that neutrinos are massive and lepton flavors are mixed. In the framework of three lepton 
families, the phenomena of flavor mixing and CP violation are described by a 3 x 3 unitary 
matrix V, which relates the mass eigenstates of three neutrinos (z^i, 2^2, 1^3) to their flavor 
eigenstates (i/g, i^/z, ^^t): 



V^i V^2 V,s \ 1^2] . (1) 




No matter whether neutrinos are Dirac or Majorana particles, leptonic CP and T violation 
in normal neutrino-neutrino or antineutrino-antincutrino oscillations depends only upon a 
single rephasing-invariant parameter J' [6], defined through 

Im(K,y^,K;.y;j = J J^ie^fs^e^j,) , (2) 

where Greek and Latin subscripts run respectively over (e, r) and (1, 2, 3). A challenging 

task of the future long-baseline neutrino oscillation experiments is to measure J', so as to 
establish the existence of CP violation in the lepton sector. In principle, a determination of 
J' is also possible from four independent moduli | ^^^1 [7], whose magnitudes can be measured 
in some CP-conserving processes. 

The unitarity of V implies that its nine matrix elements are constrained by two sets of 
normalization conditions and two sets of orthogonality relations: 

i 

E (Ki^cj) = 5,j . (3) 

a 

The six orthogonality relations define six triangles in the complex plane, as first discussed in 
Ref. [8]. Six unitarity triangles have 18 different sides and 9 different inner angles, but their 
areas are all identical to J /2. So far, some particular attention has been paid to triangles 
A3 [8] and [9] shown in Fig. 1. Because current experimental data indicate that V has 
a nearly bi-maximal mixing pattern with \Vez\ <^ 1, one can easily observe that three sides 
of A3 are comparable in magnitude; i.e., |VeiV^2l ~ I^i^fr2l ~ I^T-i^r2l ^ O-^- is therefore 
possible to establish A3 and determine its three angles [7] , once its three sides are measured 
to an acceptable degree of accuracy. In contrast, one side of A^ is much shorter than its 
other two sides; i.e., ~ 1^2^(^*21 ^ |K3V](^3|- To establish A^ needs much more 

precise data on its three sides, which must be able to show |T4iV^i| + |K:3^*3l > |K2V^2l 
l^e2V'^*2l + |K;3V^^*3l > [9]- Such au accuracy requirement is practically impossible to 

be satisfied in the near future [10]. 

To measure J and iV^o-il in realistic long-baseline experiments of neutrino oscillations, 
the terrestrial matter effects must be taken into account [11]. The probabilities of neutrino 
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oscillations in matter can be expressed in the same form as those in vacuum, however, if we 
define the effective neutrino masses and the effective lepton flavor mixing matrix V in 
which the terrestrial matter effects are already included [12]. In analogy to the definition of 
J in Eq. (2), the effective CP- violating parameter J in matter can be defined as 

ImiV^S'p.K^V;,) = J E(ea/3-y%- J , (4) 

7,fc 

where (a,/?, 7) and {i,j,k) run respectively over {e,iJ,,T) and (1,2,3). One may similarly 
define the unitarity triangles in matter with the help of the unitarity conditions 

E = Sap , 

i 

E {K^VaJ) = 5,, . (5) 
a 

For example, the effective counterparts of triangles A3 and A,- are denoted respectively 
by A3 and A,-. Because of the terrestrial matter effects, the shapes of A3 and are 
likely to be dramatically different from those of A3 and A^. It is then possible to have 
l^ei^^Til ~ I ^62^^*2! ~ iKs^/TsI some proper values of the neutrino beam energy. If this 
speculation is really true, one will be able to calculate J by using three sides of A,-. 

The purpose of this article is to carry out a systematic analysis of leptonic unitarity 
triangles in matter. We shall derive the exact analytical relations between \Vai\^ and \Vai\^ 
for a constant matter density profile. The sides of Aj (for i = 1, 2, 3) and Aq, (for a = e, yU, r) 
can then be linked to those of Aj and A^. The inner angles of Aj and A^ will also be 
calculated in terms of the inner angles of Aj and Aq,. In addition, we shall discuss the 
matter-modified rephasing invariants of V (such as J and the off-diagonal asymmetries of 
V) and illustrate the typical shape changes of Aj and Aq, with the neutrino beam energy in 
a long-baseline experiment. Our results are expected to be very useful for a complete study 
of lepton fiavor mixing and CP violation in the era of precision measurements. 

The remaining parts of this article are organized as follows. In section II, we outhne the 
master formulas to derive fki and V in terms of mj and V. Section III is devoted to the 
calculation of |V^ip. The analytical relations between the inner angles of Aj (or A^) and A, 
(or Aq,) are presented in section IV. Section V is devoted to illustrating the shape evolution 
of Aj and Aq with the matter density. Finally, a brief summary is given in section VI. 



II. FRAMEWORK 

In the fiavor basis chosen in Eq. (1), the effective Hamiltonians responsible for the 
propagation of neutrinos in vacuum and in matter can respectively be written as 
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and their difference 

n-n = 



/a 0\ 



Vo 0/ 



(7) 



signifies the matter effect [13], where E denotes the neutrino beam energy, a — y/^GpNe 
measures the charged- current contribution to the u^e" forward scattering, and A^e is the 
background density of electrons. In writing out the expression of Ti, we have assumed 
that the matter density profile is constant (namely, Ne = constant). Such an assumption is 
actually close to reality for most of the proposed terrestrial long-baseline neutrino oscillation 
experiments [14]. In order to establish the analytical relationship between V^i a-nd V^j, we 
define two quantities 

Pa/3 = '^EHafB , 

q^p = {2Ef n-UeiH (8) 
in vacuum and their effective counterparts 

~q,0 = {2Ef H-Ueifi (9) 

in matter [15], where the subscripts a and (3 run over e, and r. The determinants of Ti 
and fi in Eqs. (8) and (9) are simply det Ti = mlmlml/ {2Ey and detH = ih\fh\fh\/ {2E)^. 
In terms of (or mf) and VaiV^^ (or VaiV^j), one can obtain 



o 
1=1 

Pap = E {rf^'Va^V^^) ; (10) 

i=l 

and 

1 ^ 

^ k=l 
1 3 

= 9 E ("^i "^?^"fcV^^*fe4fe) ■ (11) 
^ k=l 

Although the exact analytical relations between Vai and Vai have been derived in Ref. [12], 

they are not simple enough to calculate VaiVp^ (or VaiV*j), which are directly relevant to the 
leptonic unitarity triangles in matter. Hence we shall establish the relations between V^jVg* 
and VcjVg* in a different and simpler way. Our strategy is as follows: first, we express VmVg* 
in terms of {Pap^'iap) ^'^'^ '^aiVpi in terms of {Vapi^ap)'^ second, we find out the relationship 
between {Pap^Qap) {PaiBTQais)'^ finally, we derive the direct relations between VaiV^^ and 
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Eqs. (3), (5), (10) and (11) allow us to express VaiV^^ in terms of {pafSilap) VaiV^^ 
in terms of {pap-i^ap)- To see this point more clearly, we introduce the coefficient matrices 
O and O : 
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(12) 



Because three neutrino masses are not exactly degenerate, the inverse matrices of O and O 
exist. Then one can obtain 







^ 5a(3 
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\Qap 



(13) 



Taking account of Eqs. (6) and (7), we find that {p^p^ Qap) (Pap^ lap) connected with 
each other through 
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(14) 



where A = 2Ea. A combination of Eqs. (13) and (14) will lead to the direct relations 
between VaiV^^ and VaiV^^, as one can sec in the subsequent sections. The usefulness of 
Eqs. (10)-(14) for discussing the matter-induced properties of lepton flavor mixing and CP 
violation has partly shown up in the literature (see, e.g., Refs. [7,15-17]) and will be further 
demonstrated in the following. 

It is worth mentioning that the generic results obtained above are only valid for neutrinos 
propagating in vacuum and in matter. As for antineutrinos, the corresponding formulas can 
straightforwardly be written out from Eqs. (6)-(14) through the replacements V =^ V* 
and A =^ -A. 



III. MODULI (a = (5) 

We first derive the relations between \Vai\^ and \Vai\'^ by taking a — (3 lor Eq. (13). 
Once Eq. (14) is taken into account, we may concretely obtain 
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where T is defined as 
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In Eq. (15), the inverse matrix of O reads 



^^2 



;^,2^ 



m2A3i(mf + m3 



^^2^ 



mfA23 

miA3i 



.12^23^31 ymiAi2(7ft? + mi) miAi2 




where Ajj = mj — rri^. With the help of Eqs. (3) and (5) as well as the relationship [7] 



Z.2 
i 



i=l 



we solve Eq. (15) and arrive at 
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A11A21 
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A12A13 


A12A13 


\Ve2\'^ 


A12A32 
A12A32 


\Ve2 




A13A22 
A21A23 






A22A32 
A21A23 


\Ve3\' = 


A23A13 
A23A13 


IK3 




A21A33 
A31A32 


iKl 




A33A13 
A31A32 



and 



K2? 



_ A31A21 




A32An 


\V,2? + 


AA13 


\Vrl 




AA23 


A31A21 


A12A13 


A12A13 


A12A13 


_ A12A32 
A12A32 


\V,2? + 


A13A22 
A21A23 




AA21 

A21A23 


\Vr2 




AA31 
A21A23 


A23A13 




A21A33 




^A32 


\Vr3 




AA12 


A23A13 


A31A32 


A31A32 


A31A32 



\Vr2? + C^ , 



^t3 



\' + C2, 
P + C3; 



and 
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= ^^IVrlf + ~^|K2P + T^I^MlP + + Cl , 

A31A21 A12A13 

~ 1 2 _ A12A32 1 2 I A13A22 , 1 2 I ^""21 1 2 I """^ai |2 ^ 

' "^^1 ~ A A ' "^^l A A ' "^^1 A I'^/^^l + \Vfj,3\ + O2 

A12A32 A21A23 

iKsP = t^^lKsr + t4^l^.iP + T^I^MaP + + ^3 , (21) 

A23A13 A31A32 

where = m? — m^, A^j = m? — m^, and 

Ci = 

C2 = 



AA13 


I^Mll 


P + 


^A23 


A12A13 


AioAiQ 


AA21 

A21A23 


IK2I 


r + 


>1A31 

A21A23 


^A32 


I^M3| 




AAr2 


A31A32 


A31A32 



A31A23 


+ A31A32 


+ A31A33 




A12A13 




A12A31 


+ A12A13 


+ A12A11 




A21A23 




A23A12 


+ A23A21 


+ A23A22 


A31A32 



^ _ ^23^12 ^ ^23^21 ^ ^23^22 

G3 . {22} 

^31^32 

In appendix A, we list the exphcit expressions of Aij and Ajj in terms of Aj^ and A. It is 
then possible to evaluate the deviation of iKip from iKip by using the formulas obtained 
above. Of course, \Vai\^ — \Vai\^ holds at the limit A ^ or £■ — > 0. 

Two sets of normalization conditions given in Eq. (5) allow us to define two off-diagonal 
asymmetries of V: 



Al = 


\Ve2\'- 




IK3P - 


IK2P = 




\Ve3\ 


Ar = 


\Ve2\'- 






\Vr2\' = 


\Vrs\'- 





(23) 



They can straightforwardly be calculated by using Eqs. (19), (20) and (21). Comparing 
between Ah (or Ar) and its counterpart Ai, (or Aji) in vacuum, we may examine the matter 
effect on the geometric structure of V. In particular, the six matter-modified unitarity 
triangles will reduce to three pairs of congruent triangles [18], if Al = or A^^ = holds. 

The sides of six unitarity triangles can also be obtained from Eqs. (19), (20) and (21). If 
three sides of a specific triangle are comparable in magnitude, one may use them to calculate 
the area of this triangle - namely, the CP- violating invariant J^. Indeed, is given by 

-iK^r - iV^ftf - |K,f - \Vpi\^y , (24) 

in which a ^ (5 running over (e,/x, r) and i ^ j running over (1,2,3). The implication of 
this result is obvious: important information about leptonic CP violation can in principle 
be extracted from the measured moduli of four independent matrix elements of V . 



IV. ANGLES (a ^ P) 

Now let us focus our attention on the inner angles of six unitarity triangles. There 
are only 9 independent angles, which appear either in triangles Ae,^,r or in triangles Ai 2,3 



7 



(see Fig. 1 for illustration). Hence it is only necessary to consider (Ag, A^, A-r) and their 
effective counterparts (Ag, A^, A^) in the calculation of 9 angles. To be specific, an inner 
angle of Aq, or Aq, (for a = e, /x or t) can be defined as 



:/3 = arg 



= arg 



Pi 



(25) 



where {a,P) run over {e,iJ,), {h,t) and (T,e), and run over (1,2), (2,3) and (3,1). 

Taking account of Eqs. (2) and (4), we obtain 



cot<l>%^jRe{Va,Vf,jV:^V;,) , 



J 



Note that J and J are related with each other through the equation [19] 

13^23 — 13^23 ■ 



(26) 



(27) 



In order to hnk cot to cot via Eq. (26), we need to find out the relationship between 
ReiVaiVpjKjVp*,) and Re{V^,Vp,V:^Vp*^ . 



For ay^(3, Eqs. (13) and (14) yield 

0-^0 



Ve2v;^ 

K2K*2 

o-'o I v^2v;. 



- AO-^T 



v,2v;^ 

^.3^3, 



(28) 



where T has been defined in Eq. (16). With the help of Eqs. (17) and (18), one may solve 
Eq. (28) and obtain the matter-induced corrections to T4jV|j. Three sides of the effective 
unitarity triangle Ag, A^ or A,- are then given by 



Ml 



A21A31 I ^11^32 .^ 



ye2v;2 = 



A21A31 

A32A21 
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A13A23 

^13^23 



Ve2V:, 



fM2 



A12A13 

A22A31 
A12A23 
A33A21 

^13^23 



Ml 



(29) 
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for A^; and 



^21^31 ^12^13 

^32^21 ^12^23 

V^aK^a = ^^yr.V:, + (30) 

^13^23 ^13^23 



for A^; and 





(A21 + ^) A31 


A21A31 




(A32 + A) A21 


A32A21 


^^3^3 = 


(Ai3 + A) A23 


A13A23 



V V* ^ (All + ^)A32 . 

^12^13 
^12^23 

v;3V;.+%±^^„,Kl (31) 

^13^23 

for Ag. It is remarkable that Eq. (29), (30) or (31), together with Eqs. (2) and (4), can be 
used to derive Eq. (27). The relationship between Re{VaiVpjV*jV/3i) and Re{VaiVpjV*jVpi) 
can also be derived from these equations. Then we are able to estabhsh the direct connection 

between cot^^g and cot0^^. 

Comparing between the definition of and the simpler notation of nine inner angles 
in Fig. 1, we have Zl = 01^, Z2 = 0^3^ ^3 = 03i^; ^4 = 0i2^ ^5 = 023^ ^ ^31.^^7 ^ ^12^ 

Z8 = 4>l^, Z9 = 0g^. One may use the similar notation (Zl, • • •, Z9) to replace (p^^^ for the 
matter-modified unitarity triangles. Therefore, 

~ (A32 + A)(A2l + ^) iAu + A}iA22+A) 

cot Zl = ^ ^ cot Zl + ^ ^ cot Z2 

A32A21 A12A12 

(A21 + A)(A22^+ A)A,3 ^3 (An + A)(A32 + ^) 2 ^ 

A12A23A12 JA31A12 ' 

~ (Ai3 + A)(A32 + ^) (A22 + A)(A33+^ 

cot Z2 = ^ cot Z2 + ^ cot Z3 

A13A32 A23A23 

A23A31 A23 ZT" A12A23 

cot 23 = + cot Z3 + (^n4-^4)(A33 + A) 

(A„ + A)(A.3 + ^)A3. (aL +^)(A3. + ^) 

Ai2A3,A3, JA23A31 I "1 ■ 

and 
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,7/1 ^32^21 A11A22 A21A22A13 A11A32 I ^|2 

cot Z4 = ^ — cot Z4 H ^ — cot Z5 H ^ — cot Z6 H ^ — \K2Ve2\ > 

A32A21 A12A12 A12A23A12 A31A12 

, 7r A13A32 A22A33 A32A33A21 A22A13 2 

cot Z5 = — cot Z5 H ^ — cot Z6 H — cot Z4 H — | VvsVggI , 

A13A32 A23A23 A23A31 A23 A12A23 

, 7c ^2lAi3 A11A33 , , . , A11A13A32 A21A33 2 /ooN 

cot Z6 = — cot Z6 H ^ — cot Z4 H ^ — cot Z5 H — VViK 1 ; (33) 

A21A13 A31A31 A12A31A31 JA23A31 

and 

,7-, A32A21 A11A22 A21A22A13 A11A32 2 

cot Z7 = ^ — cot Z7 H ^ — cot Z8 H ^ — cot Z9 H ^ — \^e2Vf^2\ > 

A32 A21 A12A12 A12 A23A12 A31 A12 

.70 ^13^32 , ,0 , ^22A33 A;-j2 A:-..--. A21 A22A13 t/* |2 

cot Z8 = — cot Z8 H — cot Z9 H ^ — cot 17 H ^ — r es Val > 

A13A32 A23A23 A23A31 A23 i7 A12A23 

,7^ A21A13 A11A33 A11A13A32 A21A33 2 /Q/|\ 

cot Z9 = — cot Z9 H — cotZ/H — cot Z8 H — Vei\/,i . (34) 

A21A13 A31A31 A12A31A31 JA23A31 

These exact analytical results clearly show how nine inner angles of six unitarity triangles 
get modified by the matter effects. 



V. ILLUSTRATION 



We proceed to numerically illustrate the matter effects on the shapes of six unitarity 
triangles, the rephasing invariant of CP violation, and the off-diagonal asymmetries of V. 
In view of current solar [1] and atmospheric [2] neutrino oscillation data, we typically take 
A21 8 X 10"^ cV^ and A32 ^ 2.3 x 10"^ cVl Wc also take 9^ ^ 33°, ^23 ~ 45°, 
6i-^ ^ 3° and 6 ~ 90° in the standard parametrization of V [7]. It is unnecessary to specify 
two Majorana-typc CP-violating phases of V in our calculations, because they play no role 
in the configuration of leptonic unitarity triangles [20]. The inputs taken above lead to 
J ^ 0.012, Ah ~ 0.147 and ^ —0.203. This means that V is asymmetric both about 
its Vei-Vfj,2-Vr3 axls and about its \43-V^2-K-i axis, and the area of each triangle is about 8 
times smaller than its maximal limit (i.e., J' = l/(6\/3) [6]). 

For a realistic long-baseline neutrino oscillation experiment, the dependence of terrestrial 
matter effects on the neutrino beam energy can approximately be written as A ^ 2.28 x 
10~^eV^£^/[GeV] [14]. This approximation is reasonably good and close to reahty, only if 
the baseline length is about 1000 km or shorter [21]. Of course, V = V holds in the limit of 
E = OT A = 0. Typically taking E = 1 CcV, 2 CcV and 3 GeV, we calculate the sides of 
six effective unitarity triangles in matter and show the changes of their shapes in Figs. 2-5. 
The corresponding results for 9 inner angles of ^e,n,T or A 1^2,3 are presented in Table 1, and 
the numerical dependence of Al and Ar on E is illustrated in Fig. 6. Some comments 
and discussions are in order. 

(a) The smallness of 613 (or {VesD Kiakes one side of A 1 or A2 strongly suppressed. When 
the terrestrial matter effect with E > 1 GeV is taken into account, three sides of Ai and 
A3 become comparable in magnitude for neutrinos {+A and V); or three sides of A2 and 
A3 become comparable in magnitude for antineutrinos {—A and V*). As a consequence of 
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matter corrections, the shortest side of A2 (i.e., S^l = jV^sKil at A = 0) turns out to be 
shorter for neutrinos; so does that of Ai (i.e., S^^ = 1^2^31 at A = 0) for antineutrinos. We 
find that the shape of A3 is relatively stable against matter corrections, no matter whether 
the neutrino beam or the antineutrino beam is concerned. 

(b) Similarly because of the smallness of ^13 (or iKsl), one side of A^ or A,- is strongly 
suppressed. The terrestrial matter effect becomes significant for > 1 GeV. In this case, 
three sides of A^ and A,- are comparable in magnitude for either neutrinos {+A and V) or 
antineutrinos {—A and V*). An interesting feature of Ae is that its side Sf^^ = \V^2Vt2\ 
dramatically sensitive to matter corrections and becomes very short for neutrinos, while its 
side Sj^]. = \VniV*i\ may significantly be suppressed by matter effects for antineutrinos. One 
can see that the shapes of A^ and A,- are relatively stable against matter corrections, no 
matter whether the neutrino beam or the antineutrino beam is concerned. 

(c) Note that the input 623 ~ 45°, which is well favored by current atmospheric neutrino 
oscillation data [2], results in \V,j,i\ ~ (for i = 1,2,3). Hence triangles A^ and A,- are 
congruent with each other; so are their effective counterparts A^ and A^-. This accidental 
result has actually shown up in Figs. 4 and 5. Prom the phenomenological point of view, it 
makes sense to measure 623 accurately and to examine its possible deviation from maximal 
mixing (i.e., 623 = 45° exactly), so as to explore the underlying /x-r flavor symmetry and its 
geometric manifestation in leptonic unitarity triangles. 

(d) Table 1 is helpful for us to understand the terrestrial matter effect on 9 inner angles 
of six unitarity triangles. One can see that Z6 and Z9, which happen to be identical as a 
consequence of 623 ~ 45°, arc relatively stable against matter corrections. In contrast, 22, 
13 and Z4 are rather sensitive to matter corrections. 

(e) No matter how the sides and angles of six unitarity triangles change with the matter 
effect, their area (J^/'i) is in most cases smaller than that in vacuum {J'/2). This unfortunate 
feature, as shown in Fig. 6, makes it hard to directly measure leptonic CP or T violation in 
any realistic long-baseline neutrino oscillations. If the neutrino beam energy is small (e.g., 
about 1 GeV or smaller), the matter-induced suppression of is not significant. In this 
case, a study of leptonic CP violation and unitarity triangles in a mec/mm-baseline neutrino 
oscillation experiment seems more feasible [22]. 

(f ) The matter effect on two off-diagonal asymmetries of V is also illustrated in Pig. 6. 
One can see that Ah ~ is a good approximation, when the beam energy of antineutrinos 
satisfies E > 1 GeV. In this interesting case, one may approximately arrive at Ag = Ai, 
A^ = A2 and A^ = A3. Such a result can also be seen from Figs. 3 and 5. The possibility 
for Ar ~ may appear only when the beam energy of neutrinos is around E 0.2 GeV, 
as shown in Fig. 6. These results are certainly dependent upon our inputs. Therefore, they 
mainly serve for illustration. 

In practice, the reconstruction of a unitarity triangle requires a series of measurements, 
which can be done in both appearance and disappearance neutrino oscillation experiments. 
More detailed analyses in this respect can be seen from Ref. [9] (see also Ref. [23]), where 
the central attention has only been paid to the triangle A,- and its effective counterpart A^. 
A similar analysis of other triangles is expected to be very lengthy and will be presented 
elsewhere [24]. As we have emphasized, it is possible to establish leptonic CP violation 
by means of the precise experimental information about the sides of six unitarity triangles. 
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In this sense, the geometric approach described here may be comphmentary to the direct 
determination of CP violation from the measurement of probabihty asymmetries between 
~^ and Va —>■ (for a ^ P running over e, /x, r) oscillations. 

VI. SUMMARY 

Considering the lepton flavor mixing matrix V and its effective counterpart V in matter, 
we have carried out a systematic analysis of their corresponding unitarity triangles in the 
complex plane. The exact analytical relations between the moduli \Vai\ and \ Vai\ have been 
derived. The sides and angles of each unitarity triangle have also been calculated by taking 
account of the terrestrial matter effect. We have illustrated the shape evolution of six 
effective triangles with the beam energy of neutrinos and antineutrinos in a realistic long- 
baseline oscillation experiment. Matter corrections to the rephasing- invariant parameter of 
CP violation and the off-diagonal asymmetries of V have been discussed too. 

We expect that this complete geometric description of matter-modified lepton fiavor 
mixing and CP violation will be very useful for the future long-baseline neutrino oscilla- 
tion experiments, although it remains unclear how far our experimentalists will go in this 
direction. We admit that how to experimentally realize the proposed ideas and methods is 
certainly a challenging question. However, it is absolutely clear that the measurement of 
leptonic CP violation must be a top task of experimental neutrino physics in the coming 
years or even decades. Let us recall what has happened in the quark sector: CP violation 
and one unitarity triangle (defined by KdKTfe + ^cdVX + VtdV^^, = in the complex plane) 
have been established at KEK-S and SLAC-S factories, and a further study of other five 
triangles will be implemented at LHC-i? in the near future. Thus we are convinced that 
similar steps will be taken for an experimental determination of leptonic CP violation and 
unitarity triangles in the era of precision neutrino physics. 

This work was supported in part by the National Nature Science Foundation of China. 
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APPENDIX A 



In the assumption of a constant earth density profile and with the help of the effective 
Hamiltonians given in Eq. (6), one may calculate the matter-corrected neutrino masses rhi 
in an analytically exact way. The relevant results have been presented in Refs. [11,12], from 
which both = mf — nij and Aij = m} — fh?^ appearing in Eqs. (19)-(22) and (29)-(34) 
can straightforwardly be read off. Explicitly, we have 

2 



A21 
A31 
A32 



Vx2-3|//3(l-;22) ^ 

2>y 2,z + ^J2,{l- z^) 



3 



(Al) 



and 



An 
A22 
A33 



1 

3 

1 

— X ^ — 
3 3 



3y 



+ ^3 (1 - 



x^ 



3t/ 



3(1-^2) 



- \zyjx'^ - 3y + A31 ; 



+ A 



21 



(A2) 



together with 



A21 = A21 + An , 
A31 = A31 + An , 
A32 = A32 + A22 , 



(A3) 



where Ajj = — rrij, and x, y and z are given by 



X = A21 + A31 + A , 
y = A21 A31 + A 



z = cos 



- arccos 



A21 (l-|K2|')+A3i (l-lKsP) 

2x^ - 9xy + 27AA21A31IT41I 



2 (x^ - 3yf^ 



(A4) 



with A = 2Ea being the matter parameter. 
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TABLES 



TABLE I. Numerical illustration of the terrestrial matter effect on nine inner angles of six 
effective unitarity triangles (Ae,//,T and Ai 2,3 in a long-baseline neutrino oscillation experiment, 
where v stands for neutrinos and V represents antineutrinos. 
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FIGURES 




FIG. 1. Leptonic unitarity triangles in the complex plane. Each triangle is named by the index 
that does not manifest in its three sides. 
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FIG. 2. The shape evolution of three effective unitarity triangles (Ai,A2,A3) with the 
beam energy of neutrinos {+A and V) in a realistic long-baseline oscillation experiment, where 
= iV^aiV^jl (for a = e,fi,T and i,j = 1,2,3) has been defined. 
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FIG. 3. The shape evolution of three effective imitarity triangles (Ai, A2, A3) with the beam 
energy of antineutrinos {—A and V*) in a realistic long-baseline oscillation experiment, where 
S^J^ = [Vai^ajl (^o^ a = e,iJ,,T and i,j = 1,2,3) has been defined. 
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V = V E=lGeY E = 2GeV E = ^GeV 
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FIG. 4. The shape evolution of three effective unitarity triangles {Ae, A^, A^) with the 
beam energy of neutrinos {+A and V) in a realistic long-baseline oscillation experiment, where 
'^ap = l^ai^^il (for ^ = 1,2,3 and a,P = e, /U,t) has been defined. 
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FIG. 5. The shape evolution of three effective unitarity triangles (Ag, A^, A^-) with the beam 
energy of antineutrinos {—A and V*) in a realistic long-baseline oscillation experiment, where 



— IT/ -V* 



(for i = 1, 2, 3 and a, fj = e, /n, r) has been defined. 
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FIG. 6. Terrestrial matter effects on J ^ A\j and for neutrinos {y with and and 
antineutrinos {p with —A and F*) in a reahstic long-basehne oscillation experiment. 
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